Tunneling dynamics of few bosons in a double well 
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We study few-boson tunneling in a one-dimensional double well. As we pass from weak interactions to the 
fermionization limit, the Rabi oscillations first give way to highly delayed pair tunneling (for medium coupling), 
whereas for very strong correlations multi-band Rabi oscillations emerge. All this is explained on the basis of 
the exact few-body spectrum and without recourse to the conventional two-mode approximation. Two-body 
correlations are found essential to the understanding of the different tunnel mechanisms. The investigation is 
complemented by discussing the effect of skewing the double well, which offers the possibility to access specific 
tunnel resonances. 
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I. INTRODUCTION 

Using ultracold atoms, it has become possible to study 
hallmark quantum effects — such as tunneling- — at an unprece- 
dented level of precision and control iH 0. S 01 • One prime 
example is tunneling of matter waves, where Bose-Einstein 
condensates have facilitated the observation of Josephson os- 
cillations Js J^|3l and the complementary nonlinear self- 
trapping IsIMIMI- In the case of Josephson oscillations, the 
atoms — initially prepared mostly in one well — simply tunnel 
back and forth between two potential wells in analogy to a 
current in a Josephson junction. However, above a critical in- 
teraction strength, the atoms essentially remain trapped in that 
well for the experimental lifetime even though they repel each 
other. 

While the above effects have been observed for macro- 
scopic coherent matter waves, many tools such as optical lat- 
tices have promoted a trend to study smaller systems with few 
atoms only. Permitting a high degree of control, they offer 
the chance to study finite-size effects and this way allow for 
a deeper understanding of the microscopic mechanisms in ul- 
tracold atoms. As an example, the recently evidenced sta- 
bility of repulsively interacting atom pairs as they move in 
a lattice lUOll . as well as the direct observation of their first- 
and second-order tunneling dynamics lITTll . should be seen as 
few-body counterparts of the above self-trapping transition. 
This motivates a thorough theoretical investigation of the few- 
boson tunneling mechanisms. 

However, while those effects are confined to the regime 
of relatively weak interactions, interatomic forces can be ad- 
justed experimentally over a wide range, e.g., by exploiting 
Feshbach resonances lfl2ll . In particular, it is well known 
that in one dimension (ID) one can tune the effective inter- 
action strength at will via a confinement-induced resonance 
lfl3ll . which makes it possible to explore the limit of strong 
correlations. If the bosons repel each other infinitely strongly, 
they can be mapped to noninteracting fermions [14], in that 



the exclusion principle serves to mimic the hard-core interac- 
tion. While the bosons share local aspects with their fermionic 
counterparts, nonlocal properties such as their momentum dis- 
tribution are very different. Sparked also by the experimental 
demonstration lUSl [la, jlTll . this fermionization has attracted 
broad interest (see lH Hi US Im IH |2l Hi IH IH and 
Refs. therein). 

In this light, the question naturally arises whether the no- 
tion of tunneling can be pushed to the strongly interacting 
fermionization limit. Indeed, a recent study has shown that 
a fermionized atom pair tunnels coherently almost like a sin- 
gle atom IztIi . In this paper, we give a systematic account of 
how few-boson tunneling evolves in the crossover from weak 
to strong correlations. Moreover, we extend that study to two- 
atom tunneling resonances occuring in asymmetric wells. 

Our paper is organized as follows. Section HI] introduces 
the model and briefly reviews the concept of fermionization. 
In Sec. [nil we give a concise presentation of the computa- 
tional method. The subsequent section is devoted to the re- 
sults on tunneling in a symmetric double well for two atoms 
(Sees. HVAWVCI i and more atoms (Sec. HVDl i. FinaUy, we 
illuminate the effect of tilting the double well in Sec. |Vl 



II. THEORETICAL BACKGROUND 
A. Model 

The subject of this article is the double-well dynamics of 
few atoms [N ~ 2 — 4), which shall be described by the many- 
body Hamiltonian (see ll24ll for details) 
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Here the double-well trap U{x) ~ + hSw{x) is mod- 
eled as a superposition of a harmonic oscillator and a central 
barrier shaped as a Gaussian S^,{x) = e-'^''^^'"^ /V2^w (of 
width w ~ 0.5, where harmonic-oscillator units are employed 
throughout.) The effective interaction in ID can be repre- 
sented as a contact potential ifisll . but is mollified here with 
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a Gaussian Scr=o.Q5 so as to alleviate the well-known numeri- 
cal difficulties caused by the S function. We focus on repulsive 
forces, i.e., g G [0, oo). 

To prepare the initial state ^'(0) with a population 
imbalance — in our case, such that almost all atoms reside in 
the right well only — we make that side energetically favorable 
by adding a linear external potential —d-x (with sufficiently 
large d ^ 0.1 — 1, depending on and g) and let the system 
relax to its ground state To study the time evolution in 

the symmetric double well (Sec.llVli. the asymmetry d will be 
ramped down to d ^ nonadiabatically (we typically choose 
a ramp time r ~ 1). By extension, it is possible to take any 
final asymmetry limf_,oo d{t) ^ 0, which allows us to look at 
the case where one well is energetically offset (SeclVTi. 



B. Fermionization 

A peculiarity of ID systems is that bosons with infinitely 
strong repulsive point interactions, g — > oo, become impene- 
trable. Mathematically, this means that its configuration space 
becomes disconnected into regions {xi ^ Xj Vi < j}, a 
feature which allows the system to be solved exactly via the 
Bose-Fenni map lfl4ll that establishes an isomorphy between 
the exact bosonic wave function ^'(J_,oo '^^at of a (spin- 
polarized) non-interacting/erm!o«;c solution ^'q , 

*i=^*o, (1) 

where A — ni;<j sgn(a;i ~Xj). The mapping rests on general 
grounds and is valid for both stationary and explicitly time- 
dependent states. Since = 1, their (diagonal) densities 
as well as their energy E will coincide with those of the cor- 
responding free fermionic states. That makes it tempting to 
think of the exclusion principle as mimicking the interaction 
{g oo), which is why this limit is commonly referred to as 
fermionization. 



III. COMPUTATIONAL METHOD 

Our goal is to investigate the few-atom quantum dynam- 
ics in the crossover to the highly correlated fermionization 
limit g ^ oo in an exact fashion. This is numerically 
challenging, and most studies on the double-well dynam- 
ics so far have relied on two-mode models ^ |29l [10, 
I3TI1 valid for sufficiently weak coupling. Here we adopt 
the Multi-Confi gura tion Time-Dependent Hartree (MCTDH) 
method ||32[ [33i l34ll . Its principal idea is to solve the time- 
dependent Schrodinger equation 

as an initial-value problem by expanding the solution in terms 
of direct (or Hartree) products $ j = ipj-^ (E) ■ ■ ■ ^ fj^., : 

(2) 

J 



The (unknown) single-particle functions (j = 1, . . . ,n) 
are in turn represented in a fixed primitive basis implemented 
on a grid. 

Note that in the above expansion not only the coefficients 
Aj but also the single particle functions ipj are time depen- 
dent. Using the Dirac-Frenkel variational principle, one can 
derive equations of motion for both Aj, tpj ll33ll . Integrating 
this differential-equation system allows us to obtain the time 
evolution of the system via (|2]i. This has the advantage that the 
basis {^j{ty\ is variationally optimal at each time t. Thus it 
can be kept relatively small, rendering the procedure very effi- 
cient. We stress that ^E* obeys bosonic permutation symmetry 
even though the direct-product basis does not; this is ensured 
by correct symmetrization of the expansion coefficients. 

Although designed for time-dependent simulations, it is 
also possible to apply this approach to stationary states. This 
is done via the so-called relaxation method Issll . The key idea 
is to propagate some wave function ^'(0) by the non-unitary 
e^^'^ {propagation in imaginary time.) As r ^ 00, this expo- 
nentially damps out any contribution but that stemming from 
the true ground state like e"'^™"^")^. In practice, one relies 
on a more sophisticated scheme termed improved relaxation 
||36[ I37I1 . which is much more robust especially for excita- 
tions. Here {'^\H\'^) is minimized with respect to both the 
coefficients Aj and the orbitals Lpj . The equations of motion 
thus obtained are then solved iteratively by first solving for 
A J W\t\\ fixed orbitals and then 'optimizing' Lpj by propagat- 
ing them in imaginary time over a short period. That cycle 
will then be repeated. 



IV. SYMMETRIC DOUBLE WELL 

Let us first focus on the tunnel dynamics in a symmetric 
well {d ~ 0). Our primary focus is on how the tunneling 
changes as we pass from single-particle — i.e., uncorrected — 
tunneling (g = 0) to tunneling in the presence of correlations 
and finally to the fermionization limit (g 00). It is natural 
to first look at the conceptually clearest situation where N = 2 
atoms initially reside in the right-hand well (Sec. lIV At . with 
an eye toward the link between tunneling times and the few- 
body spectrum (Sec. lIV BI i as well as the role of two-body 
correlations (Sec. lIV C] i. With this insight, we tackle the more 
complicated dynamics of = 3, 4, . . . atoms in Sec. IIVDI 

A. From uncorrelated to pair tunneling 

Absent any interactions, the atoms should simply Rabi- 
oscillate back and forth between both wells. This can be 
monitored by counting the percentage of atoms in the right 
well, pR{t) = (e(a;))*(i) = p{x;t)dx (p being the 
one-body density) or, correspondingly, the population imbal- 
ance S = PR — PL = 2pR — 1. Figure [TJ a) confirms that 
PR harmonically oscillates between 1 and 0. By contrast, if 
the atoms repel each other, then the tunneling process will 
be modified. For g ~ 0.2, one sees that the tunneling os- 
cillations have become a two-mode process: There is a fast 
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Figure 1 : (color online) Two-atom dynamics, (a) Population of the 
right-hand well over time, pnit), for different interaction strengths 

5 = (-), g = 0.2 (- - -), 5 = 4.7 (■ • ■), and = 25 ( ). 

(b) Snapshots of the one-body density p{x) for different times t in 
the strongly correlated case g = 25. (All quantities in harmonic- 
oscillator units throughout, see text. ) 
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Figure 2: (color online) (a) Single-particle spectrum {ta} of a dou- 
ble well with barrier height h = 8. (b) Two-particle spectrum as a 
function of the interaction strength g. Inset: Doublet formation with 
increasing g. 



B. Spectral analysis 



(small-amplitude) oscillation which modulates a much slower 
oscillation in which the atoms eventually tunnel completely 
(ppj^ « 0). In case g is increased further, we have found that 
the tunneling period becomes indeed so long that complete 
tunneling may be hard to observe. For instance, at g = 1.3 
(not displayed here) the period is as large as 2 x 10'^. What 
remains is a very fast oscillation with only a minute amplitude 
- this may be understood as the few-body analogue of quan- 
tum self-trapping, as will be discussed in Sec. 11 V B| As we 
go over to much stronger couplings (see g = 4.7), we find 
that the time evolution becomes more complex, even though 
this is barely captured in the reduced quantity [Fig. [it a)]. 
What is striking, though, is that near the fermionization limit 
(see g — 25) again a simple picture emerges: The tunneling, 
whose period roughly equals that of the Rabi oscillations, is 
superimposed by a faster, large-amplitude motion. This in- 
triguing result states that the strongly repulsive atoms coher- 
ently tunnel back and forth almost like a single particle. As 
an illustration, snapshots of the density at different t are dis- 
played in Fig.lHb): At t = 0, the fragmented pair starts out in 
the right well, and gradually tunnels to the left well until the 
fermionized pair state reemerges on the left at t w 106. 



In order to understand the oscillations, let us regard the 
evolution of the few-body spectrum {Em{g)} as g is varied 
(Fig. I2J5). In the noninteracting case, the low-lying spec- 
trum of = 2 atoms is given by distributing all atoms over 
the symmetric and antisymmetric single-particle orbital of the 
lowest doublet (illustrated in Fig.|2^). This yields the A^ + 1 
energies 

{£;™ = £o+toAW \m = 0,...,N}, 

where A^'^^ = ei — eg is the energy gap between these two 
orbitals or, in other words, the width of the lowest band. As- 
suming that for sufficiently small g still only A^ + 1 = 3 levels 
are populated in 'i'it) = J^m 6~*^'"*Cm4'm, then the imbal- 
ance ^(t) = (0(a;) — 6(— a;))^(j) (and likewise pr) can easily 
be computed to be 

d{t) = 6'-°^'^ cos(cjoiO + S^^^'^ cosiujut), (3) 

where tj,„„ = E,„, - E„ and 5^™") = 4(*„|e(a;)|*„)cmC„ 
is determined by the participating many-body eigenstates. 
Note that the term (mn) = (02) vanishes since, by antisym- 
metry, only opposite-parity states are coupled. At g = 0, due 
to the levels' equidistance, only a single mode with Rabi fre- 
quency woi = <^i2 = a'"' contributes. For very small inter- 
action energies compared to A^"', the equidistance is slightly 
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lifted, so that the Rabi oscillations are modulated by a tiny 
beat frequency ljqi — ^^12 (not shown). However, as the inter- 
action is increased further, the two upper lines £'12 virtually 
glue to one another to form a doublet, whereas the gap to Eq 
increases (Fig.|2j5, inset). 

This level adhesion, already calculated for < 5 in 
Ref. ll26ll . may be understood from a naive lowest-band two- 
mode model (see ^ for details): As g is increased, the 
on-site interaction energy eventually overwhelms the tunnel- 
ing energy A^^\ and the eigenstates evolve from number 
states \Nq^\ n[^^) in the delocalized (anti-)symmetric or- 
bitals (/ijj^o 1 in'^o superpositions of number states |iVL, iVR.) 



in the left/right-ZocaZ/ze^f orbitals (/^j^^p^ = i^o' T 
It goes without saying that any two such degenerate number 
states \v,N — v) 7^ \N — v,v) violate parity symmetry and 
only serve to form a two-dimensional energy subspace, which 
for nonzero A^"^ corresponds to the doublets in Fig. Ob). 

With these considerations on the weak-interaction behavior 
in mind, Eq. (O asserts that for times t <C T12 = 27r/cji2, 
we only see an oscillation with period Tqi ^ T12, offset by 
on a longer timescale modulates the slow tun- 
neling of period Ti2. For small initial imbalances, we have 
1 ,5(01) ^,5(12) I ^ \cq/c2\ ^ 1; SO for short times we observe 
the few-body analog of Josephson tunneling. In our case of an 
almost complete imbalance, in turn, j^^^^^ | dominates, which 
ultimately should correspond to self-trapping, viz., extremely 
long tunneling times. These considerations convey a sim- 
ple yet ab initio picture for the few-body counterpart of the 
crossover from Rabi oscillations to self-trapping. 

It is obvious that the two-frequency description above 
breaks down as the gap to higher-lying states melts (see 
Fig.[2j5), even though for two atoms no actual crossings with 
higher states occur, as opposed to iV > 3 ll26l[38ll . The conse- 
quences for the spectrum are twofold: (i) the quasi-degenerate 
doublet will break up again, and (ii) states emerging from 
higher bands will be admixed. For the imbalance dynam- 
ics, (i) implies that the "self-trapping" scenario will give way 
to much shorter tunnel periods again, while (ii) signifies a 
richer multi-band dynamics. An indication of this may be 
seen in Fig. [T] for g = 4.7, but it most clearly manifests to- 
ward fermionization, g ~2h. 

In the fermionization limit g — > 00, the system also be- 
comes integrable again via mapping ([T). As an idealiza- 
tion, assume that at t = we put two (noninteracting) 
fermions in the right-hand well, where they would occupy 

0, 1. Expressing 



the lowest two orbitals, namely V'lf ^ P 
this (fermionic) number state 'I'(O) = 



n 



through the single-particle eigenstates \n — {uaj})- via 



/3=0,1 "R 
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V2 



) leads to 



*(t = 0) 
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ao,ai6{0,l} 
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ao ' a\ I 



where li^^^ denotes occupation of the symmetric {ap 



antisymmetric (a^ 



0)or 

1) orbital in band /3. The frequencies 
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Figure 3: (color online) Top: Probability P2{t) of finding two atoms 
in the same well for g = 0, 0.2, 25. Bottom: Snapshots of two-body 
correlation function ^2(2^1, 2:2) at equilibrium points, = 0, for 

p = (t* = 44), g = 0.2 {U = 128), and 3 = 25 (t* = 53) -from 
left to right. 



En — En' contributing to ^ {t) follow in a straight- 



forward fashion: 



Ii 



AH) 



P 



'(/3) 



=0,±1 



(4) 

Moreover, let us focus on the imbalance dynamics. Since 
g{nn ) Q Qjjjy fQj. opposite-parity states n, n', the sum must 
contain only an odd number of terms. For the special case of 
two atoms, we obtain the simple result that the only partici- 
pating frequencies are A^*^^ (the lowest-band Rabi frequency, 
corresponding to the longer tunneling period) and A'^^ (the 
larger tunnel splitting of the first excited band). This links 
the strongly interacting dynamics to the noninteracting Rabi 
oscillations. 



C. Role of correlations 

In order to unveil the physical content behind the tunnel- 
ing dynamics, let us now investigate the two-body correla- 
tions. Noninteracting bosons simply tunnel independently, 
which is reflected in the two-body density ^2(2;!, X2)- As a 
consequence, if both atoms start out in one well, then in the 
equilibrium point of the oscillation it will be as likely to find 
both atoms in the same well as in opposite ones. This is il- 
lustrated in Fig. [3] which exposes snapshots p2{xi,X2; t^,) at 

the equilibrium points (where 0) and visualizes the 
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temporal evolution of the pair (or same-site) probability 

P2{t) = (e(a;i)e(x2) + e(-xi)e(-.T2))t 

= / p2{xi,X2;t)dxidx2- 

J{xiX2>0} 

As we introduce small correlations, the pair probability does 
not drop to 0.5 anymore - at g = 0.2 it notably oscillates 
about a value near 100%. This signifies that both atoms can 
essentially be found in the satne well in the course of tun- 
neling, which is apparent from the equilibrium-point image 
of p2- In plain words, they tunnel as pairs. At this point, it 
is instructive to revisit the eigenstate analysis above: While 
the g = eigenstates 2 are delocalized, at intermedi- 
ate g ~ 0.2 they have basically evolved into superpositions 
\Ni^ = 2, A^R = 0) ± |0, 2) of pair states localized in each 
well. In this light, the dynamics solely consists in shuffling 
the population back and forth between these two pair states. 

Figure [3] in hindsight also casts a light on the fast (small- 
amplitude) modulations of pr encountered in Fig. [TJa). 
namely by linking them to temporary reductions of the pair 
number p2- Thus it is fair to interpret them as attempted one- 
body tunneling. Along the lines of the spectral analysis above, 
this relates to the contribution from the ground state, in which 
the two atoms reside in opposite wells and which does tiot 
join a doublet. Since 6(a;i)6(a;2) + @{—xi)Q{—X2) is par- 
ity symmetric, only equal-parity matrix elements contribute to 
P2, which yields p2{t) w 1 - 2p("2) gj^^ (^02^/2)- 

It is clear that, as before, the time evolution becomes more 
involved as the interaction energy is raised to the fermioniza- 
tion limit (cf. g = 25). The two-body correlation pattern is 
fully fragmented not only when the pair is captured in one well 
(corresponding, e.g., to the upper right corner xi,X2 > 0), 
but also when passing through the equilibrium point t = 53. 
These contributions from higher-band excited states also re- 
flect in the evolution of P2{i), which is determined by the two 
modes uj± = A^°-' ± A*-^-*. Over time, p2 passes through 
just about any value from 1 (pair) to almost zero (complete 
isolation). In analogy to free fermions, it is again tempting 
to understand this involved pattern as two fermions tunneling 
independently with different frequencies. 



D. Many-body effects 

Although having focused on the case of = 2 atoms so 
far, the question of higher atom numbers is interesting from 
two perspectives. For one thing, at stronger interactions many 
results become explicitly A^-dependent, including distinctions 
between even/odd atom numbers ll24l l26ll . On the other hand, 
in a setup consisting of a whole array of ID traps like in iflSl 
[T^[l3l, number fluctuations may automatically admix states 
with A^ > 2. 



1. Complete imbalance 

For A^ > 3, the weak-interaction behavior does not differ 
conceptually. In fact, Eq. (O carries over, 

E '5''™'cos(t^„„t), 

m < n 

but with the sum now running over < m < n < A^. Strictly 
speaking, the dynamics is thus no longer determined by two 
but rather in principle N{N + l)/2 modes {ran) - although 
about half of these fail to contribute by symmetry. Nonethe- 
less, the basic pattern can be understood from the two-atom 
case, as will become clear in a moment. 

For g = 0, assume an ideal initial state with all atoms in 
the right-localized orbital 0r = (0o + 0i) of the lowest 
band. The weight coefficients cn{No) = (A^o, - ^o|*(0)> 
with respect to the eigenstates |Ao, A^i) have a binomial dis- 
tribution 

which for larger A^ asymptotically equals a Gaussian, with a 
sharp peak (AA^o = VN /2) near No = A^/2. In fliis fight, 
only these few states should contribute. Again, the equidis- 
tance of the levels guarantees a simple imbalance oscillation 
with A(°^ For interaction energies small compared to A^^\ 
the Rabi oscillations will again be modulated by beats, similar 
to the case N ^ 2. 

As we move to larger values g ^ 0.2, the higher-lying of 
the N + 1 levels have again merged into doublets luql . In 
particular, the highest eigenstate pair was conjectured to be 
roughly of the form | A^l = A^, A^r = 0) ± |0, A^) (in the fimit 
h 00). The idealized state distribution should be peaked 
at just these two vectors, whose energy splitting in the bare 
two-mode model has been estimated as w ^ 2NU / {N — 
1)! X (2A(°Vt^)^ iH, where U denotes the on-site inter- 
action energy. Thus the tunnel period is expected to grow ex- 
ponentially as A^ — > 00, a trend which may be roughly extrap- 
olated from Fig. |4] (insets). Ultimately, this should connect 
to the condensate dynamics vafid for A^ > 1 ||6l[28il29i[30ll. 
when tunneling becomes inaccessible for all intents and pur- 
poses. Of course, realistically, neighboring states will also 
be excited, which makes the time evolution richer However, 
the separation of time scales leads to the characteristic inter- 
play of fast, small-amplitude oscillations (related to attempted 
single-particle tunneling) and a much slower tunnel motion, as 
observed in Fig.|4] 

Things become more intricate if we leave the two-mode 
regime, cf. g = 4.7. As has been demonstrated in ||26ll , (anti- 
)crossings with higher-lying states (which connect to higher- 
band states at g = 0) occur for A^ > 3. Given our experience 
of the two-atom case, one might again expect a simplified be- 
havior as we approach the fermionization limit. However, we 
will argue below that this has to be taken with a grain of salt 
because an initial state with A^ hard-core bosons in one well is 
highly excited. In the spirit of the Bose-Fermi map, an ideal- 
ized state with A^ fermions prepared in one well will have con- 
tributions from all excitations llio"*; liV; ■ • ■ ! 1on_i^^)- (a^g = 
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Figure 4: (color online) Time evolution PK(t) of (a) N = 3, (b) N = 4 atoms initially in one well. Shown are the coupling strengths g — 

( — ), g — 0.2 ( ), g = 4.7 (■ • ■), and g = 25 (— • — ). Insets: Long-time behavior for g = 0.2 (the longer period) and g — 4.7. (Observe 

the different time scales in both insets.) 



0, 1 V/3) in the N lowest bands, which is proven by induction 
on = 2. In view of (|4|i, many more frequencies are ex- 
pected to be present: Besides the individual tunnel splittings 
for each band, these should in principle be all four com- 
binations A(°) ± A(i) ± A(2) for iV = 3, and 4 X 4 combina- 
tions {A^ ±AW±A(") I < / < TO < n < iV}foriV = 4 
etc, taking into account parity-selection rules. However, in 
the fermionization limit with the idealized initial state above. 
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-the 



things simplify even further. Since TVpi = 
Fock-space representation of Q{x) in the context of Eq. (|3]l- 
is a one-particle operator, an eigenstate |n)_ is coupled only 
to "singly excited" states of the type |n')_ = a^' 
(for some f3), with an excitation frequency ojr, 
This yields an imbalance of 



n. 



A(/3). 



1 



JV-l 



/3=0 



This simple formula should be contrasted with the surprising 
complexity of the fermionization dynamics already for atom 
numbers as small as = 3,4. This is illustrated in Fig. |4] 
where p^it) = [{6{t) + l]/2 is plotted (cf. g = 25). To 
be sure, for finite g and using a realistic loading scheme, a 
few more modes contribute, thus naturally rendering the dy- 
namics more irregular. But even the inocuous formula above 
can account for the seemingly erratic patterns in Fig. |4l The 
key to see this is to consider the distribution of frequen- 
cies {A^'')}. In the unreaUstic limit that A(^) « A^") V/3, 
the imbalance would be a neat Rabi oscillation for any A^, 
6{t) w cosA'^^t. However, a realistic barrier likely has a 
Gaussian-type shape and a finite height; hence the splittings of 
higher bands tend to grow monotonically. As a consequence, 
only the lower-band frequencies A^^^ will contribute to the 
tunneling, whereas the higher-band splittings make for much 
faster modulations, which average out on a larger time scale. 
The gist is that for N ^ 1, those few lowest-band modes 
only have a weight of 0(1/A^), so in a realistic scenario one 
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Figure 5: (color online) Partial-imbalance effects in the fermioniza- 
tion limit (g = 25). (a) Small-imbalance oscillations (Scenario 1.) 
for Ai" = 3, 4 atoms. Plotted is the population of the right-hand well, 
PB.{t)- Bottom: Density evolution p{x;t) for A*' — 1 = 2 (b) and 
N — 1 — 3 atoms (c) initially in the right-hand well if exactly one 
atom is present on the left (Scenario 2.). 



expects quasi-equilibration around PR = 1/2. 



2. Partial imbalance 

While we have so far assumed that all atoms are prepared 
in one well, it is natural to ask what the effect of incomplete 
imbalances pyi{Q) < 1 would be. For simplicity, we will focus 
on the fermionization limit (here g = 25). Two scenarios are 
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conceivable, in principle: 

1. Small imbalances pfj w 1/2, i.e., small perturbations of 
the ground state; 

2. Preparing, say, iV — 1 atoms in one well and one in the 
other. 

Option (1.) is plotted in Fig. Ha) for TV = 3,4. 
We clearly observe Josephson-type oscillations in each 
case, but with markedly different time scales. This may 
be understood from the spectral structure near fermion- 
ization i26ll : For even N, the fermionic ground state 
. . . , has all bands filled, so 

that the lowest excitation is created by moving one atom from 
band f3 = N/2 - 1 to /3 = N/2. Thus the "Josephson" fre- 

{N/2) {N/2-1) . , ■ , , J 

quency luqi = Eq ^ is a large mter-band gap, 

which for = 4 gives a period of Tqi w 4. For odd N, by 
contrast, the mechanism is a different one: Here the ground 
state leaves the highest band only singly occupied, so that the 
lowest excitation frequency is the small intra-band splitting 
woi = A(^-i)/2 j„ Fig.Ba) {N ^ 3), this may be identified 
as the rather long period Tqi w 40. 

Scenario (2.), paraphrased in the case = 3, is the ques- 
tion of the fate of an atom pair if the target site (the left well) 
is already occupied by an atom. The striking answer, as evi- 
denced in Fig.|5jb), is that the process can be viewed as single- 
atom tunneling on the background of the symmetric two-atom 
ground state. The tunneling frequency in the fermionization 
limit is A^^^ « 27r/40, which has the intuitive interpretation 
of a fermion which — lifted to the band /3 = 1 — tunnels inde- 
pendently of the two lowest-band fermions. From that point 
of view, it should come as no surprise that adding another par- 
ticle destroys that simple picture. In fact. Fig. [Sjc) reveals 
that if we start with — 1 = 3 atoms on the right, then the 
tunneling oscillations appear erratic at first glance, and a con- 
figuration with three atoms per site becomes an elusive event 
(see, e.g., t « 22, 44 or 72). In the fermionic picture, this can 
be roughly understood as superimposed tunneling of one atom 
in the first excited band (A^^^) and another in the second band 
(A(2) ^ 27r/15), while the remaining zeroth-band fermions 
remain inactive. 



V. ASYMMETRIC DOUBLE WELL 

We have so far used the tilt d of the double well merely as 
a tool to load the atoms into one well. The question naturally 
arises whether the actual tunnel oscillations can be studied in 
asymmetric wells so as to manipulate the nature of the tun- 
neling. Specifically, we consider a setup similar to Sec. |IV] 
Two atoms are prepared in the right well (i.e., in ground state 
^g'^°-'with a large initial asymmetry do). Subsequently, the 
asymmetry is ramped down to a final value d ^ 0, thus trig- 
gering the tunnel dynamics. 



A. Tuning tunneling resonances 

In symmetric wells, pair tunneling is always resonant in the 
sense that an initial state with all atoms on one site is equal in 
energy to one with all atoms in the opposite well iflll [ssll . 
Conversely, single-atom tunneling should only be likely so 
long as the repulsive interaction does not shift the pair state's 
energy off resonance with a target state of only a single atom 
on the left. This squares with our finding that the pair prob- 
ability p2 (Fig.[3]l drops to 50% in the equilibrium points for 
g = 0, while in the correlated case (g = 0.2) it does not vary 
considerably from unity. To condense this insight into a single 
quantity, let us define 

pi = max{l -P2(i)} 

as the (maximum) single-atom probability, relating to the 
event of finding the atoms in different wells. 

Figure |6] shows how pi changes when the final asymme- 
try d between the wells is varied. For g = 0, pi{d) has a 
plateau for d < 0.011. This relates to the transition from co- 
existence of single-atom and pair tunneling (at c? = 0) to the 
point where the right-hand well is lowered such in energy that 
the initial pair state energetically matches a state with exactly 
one atom on the left. From the perspective of the two-body 
density in Fig. [3] the final state at d = 0.011 corresponds to 
the equilibrium-point snapshot for d = 0. For larger values 
of d, the energy difference between both wells is too large to 
transfer a substantial fraction of the population to the other 
well. 

By contrast, at g ~ 0.2 the repulsion is sufficiently strong 
to drive the single-atom tunneling off resonance at d — 
(Fig.|6]l. Lowering the right well so as to compensate for the 
interaction-energy shift leads to a dramatic increase of the tun- 
nel amplitude near d = 0.038. The value of sa 1 confirms 
that this is pure single-atom tunneling: After half a tunnel pe- 
riod, both atoms are found precisely in opposite wells, until 
they return to the pair state on the right site. 

Despite the more convolved dynamics that emerges as we 
go higher interactions, the one-atom tunnel resonance persists. 
However, in the fermionization limit g oo, yet another 
resonance emerges at d = already (Fig. |6]l. As in the un- 
corrected case, this signifies coincident single-atom and pair 
tunneling. This resonance, however, is much more sensitive 
to symmetry breaking, which is intelligible from the picture 
of two fermions hopping simultaneously in different bands 
/3 = 0, 1. Skewing the double well (d > 0) thus attenuates 
both one- and two-atom tunneling until another, pure single- 
atom resonance is hit at d ~ 0.58. Conversely, energetically 
lifting the right-hand well (d ~ —0.5) makes tunneling to ex- 
cited target states accessible. 

B. Spectral analysis 

To better understand the dependence of the tunnel dynamics 
on the tilt d, let us consider the two-body spectrum {Em{d)} 
at fixed coupling g. Since both the noninteracting and the 
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Figure 6: (color online) Maximum single-atom probability pi as a 
function of the tilt parameter d. Solid line: uncorrelated tunneling, 
g = 0; dashed line: con'elated tunneling, g = 0.2. Inset: Near 
the fermionization limit, g — 25. Note that the resonances are not 
symmetric in d owing to the unsymmetric initial state ^'(0). 




0.1 0.2 0.3 0.4 0.5 0.6 
asymmetry d 



Figure 7: Single-particle spectrum {ea{d)} of a double well as a 
function of the asymmetry d. 

fermionization limit can be deferred from the single-particle 
spectrum, we will first stop to review the tilted double well. 

1. One-body spectrum 

Figure |7] displays the spectrum {ea{d)} of the double well 
U{x) = + h6yj{x) — d ■ X for variable asymmetries d. 
For simplicity, let us resort to a simple model and expand the 
one-body Hamiltonian h{p, x) ~ ip^ + U (x) in terms of two 
modes (j)a=L{R) localized on the left (right) site (tacitly as- 
suming a fixed band f3). We denote by 

• {4'a\h\(pa) = e± <;/2 the energies pertaining to /io/flfec/ 
wells, where the left site has an energy offset 

• |((^l|^|'/'r)| = a/ 2 the tunnel coupling. 

Then a straightforward diagonaUzation yields 

(Pa,, oc A • + k ± A(<^)] (a = 0,1) 

e,,, = eT^A(0 

where A(<;) = vA^-T^ is the energy gap in the presence 
of the tilt. In the symmetric case, the states are simply given 
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Figure 8: (color online) Two-body spectrum {Em,{d)} in a tilted 
double well, (a) g = 0, 0.2 (b) g = 25. 

by the (anti-)symmetric orbitals (/)a,?=o oc (pi, ± (pn, with the 
usual tunnel splitting A(0) = A. As we switch on a tilt <; > 
0, parity is broken and the once delocalized states break up 
into one decentered on the left {(pi) and one on the right {(pg) 
as <; ^ A. This goes along with a level repulsion of eo/i,? 
about <j = 0, where the 0i state pinpointed on the left site is 
energetically lifted, and vice versa. As the states decouple for 
<; ^ A, the energy approaches that of the isolated subsystem 

The above picture holds for each band /3 individually, pro- 
vided their levels are well separated. In fact. Fig. |7] confirms 
that scenario for tilts small compared to the interband gap, 

<C e'''"''^^ — i'^^^ ^ 2. For strong enough asymmetries d, 
though, states emerging from different bands mix, and new 
avoided crossings are observed in the plot. 



2. Two-body spectrum 

Noninteracting limit: In the uncorrelated system, g ~ 0, 
the many-body spectrum {En = "a^a} is obtained from 
the number states \n) of the single -particle eigenstates (pa- 
The energy shift of the levels En{d) with respect to d = 
thus depends on the balance between contributions from sym- 
metric orbitals ^q'^'' and antisymmetric ones. Specifically, the 
d = Q ground state exhibited in Fig.[8ja) is a coherently sym- 
metric state 1 20) = [00°'' ]^^- Consistently, for perturbations 
d > it localizes on the right, with its level shifting downward 
- in stark contrast to the second excitation |02) = [0^"']®^. In 
between, |11) is a compromise between these two borderline 
cases in that both partial energy shifts cancel out, leaving a 
delocalized state. This gives us a new perspective on the tun- 
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neling dynamics reflected in Fig. |6] Imagine we start with 
afl atoms prepared in the right wefl, viz., the ground state 
^(d^oc)^ and then ramp down d{t) so as to trigger 
the tunnehng. If we follow the ground-state level nonadia- 
batically, then at d = it finds three closely packed levels 
EjniO) it can couple to - in the sense that | {"^^^1^(0)) I 7^ 0, 
so that a nontrivial dynamics becomes possible. In fact, at 
d = 0, these correspond to Rabi oscillations. If we were to 
choose a final asymmetry d ^ 0.01 (in the notation above, 
<^(o) < a'")), roughly the same level would be available, con- 
firming the plateau encountered in Fig.|6] However, for final 
values d > 0.01, the levels decouple, and no longer are there 
any target states at disposal for tunneling. 

Medium interactions: These elementary thoughts also 
help us explore the nontrivial dynamics for intermediate cou- 
plings, as shown for g = 0.2 in Fig. |8ja). The d = 
ground state, in the limit A'"' 0, has the Mott-insulator 
form |1l1r) and should be insensitive to symmetry break- 
ing d > 0. By contrast, the quasi-degenerate excited pair 
|2l0r) ± |0l2r) only requires a minute perturbation to break 
up into two localized states. It is plain to see that, at d « 0.04, 
the lower excited curve anti-crosses the ground state, and the 
two states are virtually swapped. Resorting again to a simple 
two-mode model, the (avoided) crossing occurs for tilts (, ~U 
matching the on-site repulsion energy. 

The bearing this has on the tunnel dynamics is evident: 
Apart from the self-trapping scenario at d = 0, there is a 
fairly broad tunnel resonance at d « 0.04, where the fully 
imbalanced initial state ^'(0) couples to that with one atom 
on each site, |1l1r)- This is but the one-body resonance en- 
countered in Fig. |6] To come by a crude estimate for the crit- 
ical value dc, assume that the energy of initial and final states 
match, {Hd^)i — {Hdjt- Modeling the initial pair state by the 
ground state ^'q''"'' (at the initial do > 0), and the final state 
with a single atom on the left by ^'q^'', yields the estimate 

in terms of the ground-state energies at the initial do > and 
d = 0, respectively, and the elongation (x) at time t = 0. 

Fermionization limit: Figure |8jb) shows the spectrum 
near fermionization, g = 25. The d = ground state turns out 
to be widely robust against perturbations, which can be under- 
stood from the fact that its fermionic counterpart |1q°-'iJ*''')_ 
has balanced populations of right- and left-localizing orbitals. 
The only way to obtain a right-localized ground state is to 
lower one well enough for it to hit a localized state from the 
upper band. This is what happens at d w 0.6, where the tilt 
energy <,^^^ /2 = e'^^ — e'"' compensates the inter-band gap. 
That crossing marks just the one-body resonance seen in Fig.|6] 
at d w 0.6. In the fermionic picture invoked above, it may be 
thought of as one excited fermion tunnehng to the lowest level 
on the left. 

If we follow the localized state nonadiabatically, then at 
d = we recover the mixed single-atom/pair resonance laid 
bare in Fig. |6] Further ramping up the right well to d w —0.3 
(where the spectrum is mirrored at d = 0), we see yet an- 
other crossing. A closer look reveals that the partner state is 



entirely localized on the left, so that one might hope for a pair 
resonance. However, as both states are localized in disjoint 
regions, they are not coupled by the perturbation (— d • x), and 
in practice no tunnel resonance is observed. It may be illu- 
minating to look at this from the fermionic perspective. For 
d « -0.3, the initial state on flie right is ^-(0) « |1^°); lj^^)_, 
while the partner state emanating from £'(0) « 8 in turn is 

given by lip"''; Ip^'')-- In this light, the tunneling "resonance" 
in question refers to the following situation: Two fermions si- 
multaneously hop from the zeroth (first excited) level on the 
right down to the zeroth level (up into the second level) of 
the energetically lower left site. While both processes indi- 
vidually are off resonance, the total energy is conserved. This 
reflects in the one-body spectrum (Fig. where no avoided 
crossing is to be observed at d « —0.3 - rather, there is an 
accidental crossing of the sums En = "a^a- However, 
at d w —0.6, another avoided crossing emerges, which — in 
the fermion language — corresponds to multiple one-body res- 
onances with the first and second excited level in the left well. 



VI. CONCLUSIONS AND OUTLOOK 

We have analyzed the crossover from uncorrected to 
fermionized tunneling of few ID bosons in a double well. 
The pathway leads via strongly delayed pair tunneling for 
medium interactions — associated with doublet formation in 
the few-body spectrum — to fermionized tunneling, where the 
strongly correlated atoms tunnel back and forth with charac- 
teristic modulations. By analogy to free fermions, these may 
be understood as multi-band Rabi oscillations, which become 
more and more complex and quasi-equilibrate for large atom 
numbers. To uncover the physical mechanisms, it is essential 
to study two-body correlations. These reveal a strong suppres- 
sion of single-atom tunneling for intermediate coupling, with 
a revival toward fermionization, where an involved interplay 
of pair and single-atom tunneling is observed. 

Whereas for small interactions, higher atom numbers es- 
sentially only increase the tunnel period but do not change 
the scenario qualitatively, the multi-atom dynamics becomes 
much richer as fermionization is approached. Apart from the 
above case of a complete initial imbalance, this applies to situ- 
ations where not all atoms are initially in one well. In particu- 
lar, Josephson-type small-amplitude oscillations exhibit vastly 
different time scales for odd/even numbers. On the other hand, 
initially storing an extra atom in the target well suppresses the 
lowet-band tunneling and thus leads to a simplified dynamics. 

Finally, studying the dynamics in asymmetric wells pro- 
vides a valuable perspective on the tunnel mechanism in terms 
of one- and two-atom tunnel resonances. Depending on the 
energy difference between the sites, the tunnel amplitude can 
be largely enhanced or suppressed. For noninteracting bosons, 
this has been described by a plateau of the single-atom prob- 
ability about the asymmetry parameter d = 0. At medium 
interactions, in turn, single-particle tunneling becomes reso- 
nant only when the energy offset of one well compensates the 
interaction-energy shift at d > 0. In the fermionization limit, 
another d = resonance emerges, accompanied by higher- 
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level resonances at d 7^ 0. Those features are explained in 
terms of avoided crossings in the spectrum as d is varied. 
Such a deeper understanding of the tunneling may pave the 
way to an active control of strongly correlated systems, for 
instance by allowing to transport definite numbers of atoms 
from a reservoir to a target well. 
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